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Abstract 

An abstract DNA-type system is defined by a set of nonlinear kinetic equations with polynomial 
nonlinearities that admit soliton solutions associated with helical geometry. The set of equations 
allows for two different Lax representations: A von Neumann form and a Darboux-covariant Lax 
pair. We explain why non-Kolmogorovian probability models occurring in soliton kinetics are 
naturally associated with chemical reactions. The most general known characterization of soliton 
kinetic equations is given and a class of explicit soliton solutions is discussed. Switching between 
open and closed states is a generic behaviour of the helices. The effect does not crucially depend 
on the order of nonlinearity (i.e. types of reactions), a fact that may explain why simplified models 
possess properties occuring in realistic systems. We explain also why fluctuations based on Darboux 
transformations will not destroy the dynamics but only switch between a finite number of helical 
structures. 

PACS numbers: 87.23. Cc, 03.67.Mn, 05.45.Yv 
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I. INTRODUCTION 



There are certain aspects of the DNA dynamics that go far beyond standard fields of in- 
terest of genetic engineering. A striking example is the observation that DNA is essentially a 
universal Turing machine plQ]. The hardware of the machine consists of two strands (input- 
output tapes) and an enzyme plays a role of finite control. As opposed to our PC's that 
cannot work without electricity, the DNA machine is powered by networks of autocatalytic 
reactions. The double-helix structure is a trademark of the machine. 

A new DNA helix is formed by a polymerase that moves along an existing strand and 
gradually, step by step, produces the complementary strand. This seems to be a continuous 
process that resembles the motion of finite control in a real-life Turing machine. The discrete 
time-steps of an abstract Turing machine can be obtained if one looks at snapshots of the 
evolution taken at times t n = nAt, say. Finally, the two helices start to unwind, and the 
two strands behave as the tapes of two new independent Turing machines. This process is 
called replication. 

The problem how to theoretically model the DNA dynamics is an old one. Between 
various notable theoretical constructions one should mention classical Hamiltonian systems 
(typically of a lattice type) ODD or replicator networks ODD 

The above approaches are well motivated. There is little doubt that DNA is a nonlinear 
dynamical system. It is also clear that the DNA evolution is a sequence of chemical reactions 
and we know how to translate kinetic schemes into differential equations with polynomial 
nonlinearities. Therefore no matter what networks of reactions one will invent, the result is 
in general a set of nonlinear equations for certain nonnegative variables. As opposed to the 
mechanical models, the kinetic equations are not directly linked to soliton systems. 

The DNA dynamics does not look chaotic. One could even argue that organisms whose 
DNA evolves in a chaotic way would be eliminated by natural selection. Systems that are 
nonlinear and non-chaotic are almost certainly solitonic. Solitons are localized, propagate 
in a stable way, and are very resistant to external noise. Soliton equations naturally lead 
to 'switches', that is, solutions of a kink type. This is why, the line of research that deals 
with Hamiltonian systems is concentrated on 'DNA solitons'. Another argument in favour 



of DNA solitons is that soliton scattering can serve as a model of computing 
The DNA computer might be a soliton machine. 
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The latter observation points into another general class of dynamical systems: Those 
whose discretized dynamics (in the sense of t n = nAt snapshots) is equivalent to some 
algorithm. The first example of such a link is the Toda lattice interpreted in terms of 
spectrum generating algorithms Q. Later on, the construct™ was generalized to a general 
class of 'Lax equations' |14U15I]. The term 'Lax equation' is often applied to equations whose 
right-hand side is given by a commutator. However, in the soliton case we need more than 



just a commutator: We need a Darboux covariant Lax pair ^| and this is not exactly 
synonymous to a commutator. Equations with commutator right-hand sides may not be 
Darboux covariant. Or, even if they are Darboux covariant, the Darboux transformation 
can be inconsistent with some important constraints. 

To sum up, it looks like it would be ideal to work on DNA dynamics with soliton kinetic 
equations that have some sort of Darboux-covariant Lax representation. Such equations do, 
in fact, exist and can be characterized in quite general terms. To find them we are led by 
still another hint: The ubiquity of helical structures associated with DNA dynamics. The 
fact that the helices are formed (and deformed) as a result of the dynamics suggests their 
3D parametrization of the form (t, x(t),y(t)), where t is time (i.e. position of the read/write 
head) and x(t), y(t) are related to kinetic variables. The kinetic variables must therefore 
occur in pairs. Alternatively, having an even number of dynamical variables one can always 
work with their complex representation £(t) = x(t) + iy(t). 

Let us now try to derive a simple kinetic scheme that leads to a helix. At this stage we 
are satisfied by the 'boring' case of a helix that does not replicate. Consider the following 
catalytic network [12 1 

B ^ A + B (1) 
A + X X+ (products) (2) 
A + B A + (products) (3) 

Y ^ B (4) 

Assume the reactions (j2J, (jSJ) are 0-th order in A and B, respectively. Then 

[A] = k[B] - k>[X\, [B]=k"[Y]-k[A], [X] = [Y] = 0. (5) 

If we define two new (in general non-positive) variables x — [A] — (k"/k)[Y\, y = [B] — 
(k'/k)[X], we can rewrite the system as i£ = k£ where £ = x + iy. The system evolves as a 



harmonic oscillator. The curve 1 i— » (t, [A(t)], [B(t)}) is a helix. If 

-(::)■'-(;:)■ 

a>0,b>0,a + b>0,a = b = 0, then ip = [H, p\. This is what we call a Lax form of the 
kinetic equations. If eigenvalues of p are nonnegative at some t = t Q then for any projector Pa 

(A (i 

we have [A] = Tr pP^ > 0. Let now Pa and Pb project, respectively, on and 

W V- 

Thenx= [A]-(a+b)/2 and y = [B]-(a+b)/2, i.e. (ife'/Jfe)^] = (Jfe"/A;)[y] = (a+6)/2. The 
kinetic variables satisfy = Tr p(t) Pa and [-B(t)] = Tr p(i)Pg. The condition [.£?(£)] > 

can be always satisfied by an appropriate shift of origin of the (x, y)-plane. What is 
interesting, the projectors Pa and Pb do not commute (projectors in a plane commute if 
and only if they project on either parallel or perpendicular directions, which is not the case 
here) . This type of noncommutativity is an indication of a non-Kolmogorovian probabilistic 
structure behind the kinetics. 

The above analysis explains the motivation for the formalism we develop in this paper. 
We will work with kinetic equations that possess a two-fold Lax structure. First of all, all 
these equations can be written as von Neumann type nonlinear systems. Secondly, each of 
these equations can be regarded as a compatibility condition for a Darboux-covariant Lax 
pair. 

The Darboux transformation will also play a two-fold role in our formalism. It will allow 
us to find nontrivial exact solutions of systems of coupled nonlinear kinetic equations. The 
transformation is known to switch between different solutions of the equation in a way that 
respects certain conservation laws (Casimir invariants). It is therefore also a natural formal 
object representing fluctuations along a symplectic leaf of the dynamics. This effect will 
lead to both mutation and formation of an open state in the associated helical structures. 

That nonlinear von Neumann equations form a dynamical framework for various tradi- 
tionally non-quantum domains was the main message we tried to convey in In the 
present paper we lay special emphasis on the aspect that was implicitly present in all the 
examples we gave, but which at a first glance may not be completely obvious: The fact 
that von Neumann equations are not analogous but, in fact, exactly equivalent to systems of 
kinetic equations. Any von Neumann equation is just a form of a Lax representation of a set 
of kinetic equations (this observation was implicitly present also in the analysis of replicator 



equations given in [19j . although technically this argument was not equivalent to the one we 
give below). This is why any integrable dynamical system that involves kinetic variables is 
suspected of being implicitly von Neumannian. 

In spite of its simplicity, the above statement is surprisingly difficult to digest by the 
community of chemists or molecular biologists. The roots of the difficulty seem to be related 
to the misleading term 'quantum probability' suggesting that non-commutative propositions 
are restricted only to microscopic systems. 

Therefore, one element we will need to understand before launching on technicalities 
of soliton kinetic equations, is the structure of probability models that are behind soliton 
kinetic evolutions. We have seen already at the simple example discussed above, that the 
'concentrations' [A] and [B] were linked to noncommuting projectors Pa and Pb- This type of 
noncommutativity is a typical feature of soliton kinetics and indicates that the probability 
calculus is here nonclassical. We therefore need to clarify some formal aspects of non- 
Kolmogorovian probability, and explain why this is the type of formalism that necessarily 
occurs in chemical kinetics. 

II. HILBERTIAN MODEL OF PROBABILITY 

A proposition is represented by a self-adjoint operator with eigenvalues and 1, i.e. a 
projector P = P^ = P 2 acting in a Hilbert space TC. Probability is computed by means of 
the formula 



where p is a 'density operator' ('density matrix') acting also in 7i and satisfying p = p^ 
(since then probabilities are real: Tr Pp = Tr Pp = Tr Pp^), p > (i.e. the eigenvalues of 
p are non-negative) since then p > 0. The normalization of probability leads to Trp = 1: If 
the set {Pj} of orthogonal projectors is complete in the sense that ^ ■ Pj = 1 then 



In a given Hilbert space there exists an infinite number of such complete sets. If a given prob- 
lem can be represented entirely in terms of probabilities associated with a single complete 
set, then all the probabilities will satisfy criteria typical of Kolmogorovian models. 



p = Tr Pp 



(7) 




(8) 
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If the kinetics is implicitly Hilbertian then all the dynamical variables, say pi, P2, ■ ■ ■ , Pn, 
have to allow for a representation of the form (JJJ) . If a given set of kinetic equations couples 
probabilities Pi,...,Pn that do not belong to the same complete set then ^2j =1 Pj 7^ 1- This 
happens whenever the set {Pj} involves projectors that do not commute with one another 
since then necessarily . Pj ^ 1. Actually, this is what we found in the simple kinetic model 
of harmonic oscillators discussed in the Introduction, where Pa and Pb did not commute 
and, hence, did not belong to the same maximal set. 

It follows that the conditions we impose on the kinetic equations can be reduced to the 
requirements imposed on the dynamics of the density operator p(t): For all times t± < t < ti 
in the interval of interest we must find (i) p(t) = p(ty, (ii) p(t) > 0, and Tr p(t) = const < 00. 

A basis in the Hilbert space is given by a set {|n), (n\m) = 5 nm }. Matrix elements of the 
density operator are in general complex 

Pnm = (n\p\m) = x nm + iy nm (9) 

Pnm •^nm ^Vnm •^mri ^Vmn Pran (10) 
and thus X mn X nm) t/mn Vnmi Vnn 0' 

The diagonal elements p nn = x nn are themselves probabilities since 

Pnn = (n\p\n) = Tr P n p = p n = x nn (11) 
where P n = \n){n\. Now, let 

i#> = ^(Ij> + i*», m = ^{\j)-m), (12) 

P jk = \jk){jk\,P; k =\jk'){jk'\. Then 

x jk = Pjk ~ ^Pj ~ ^Pk, (13) 
Vjk = P'jk ~ \P] ~ \pk, (14) 

e i7r/4 

Pjk = Pjk + ip'jk ~ ~Jf(Pj + Pk), (15) 

where pjk = Tr PjkP, p'j k = Tr Pj k p. It follows that a single density operator p encodes three 
families of probabilities: {p n }, {Pjk}, and {p'j k }- They are associated with three families of 
projectors: {P n }, {Pjk}, and {Pjk\- Additional relations between the probabilities follow 
from 

x 3i = Pi = Pn ~ I'r Vjj = = P'jj ~ Pi (16) 



and the resulting formula pj = p'^ = Pjj/2. Let us note that {P n } is complete, i.e. 
^dim/Wp^ _ ^ Yln^i^ Pn — 1- I n order to understand completness properties of {Pjk} 
and {-Pjfc} we introduce, for j < k, two additional types of vectors and their associated 
projectors: 

\jk ± ) = ^(\j)-\k)), \jk' ± ) = ^=(\j)+i\k)), (17) 

Pjk = \jk )0'& J i> = \jk' L ){ik' L \. The completeness relations for {Pjk} and {Pjk} 
follow from the formula 

p jk + p 3 i = p j + p k = p; k + p£. (is) 

The non-Kolmogorovity of the set of probabilities can be illustrated by means of uncer- 
tainty relations. Indeed, for any three operators satisfying [A, B] = iC one can prove the 
uncertainty relation AAAB > \\{C)\ where AA = a/ (A 2 ) - (A) 2 , (A) = Ti(Ap), etc. 
If A = P = P 2 is a one-dimensional projector then (P) = p is a probability and one 
finds AP = v/p(l — p). We say that two propositions Pi and P2 are complementary if 
AP l AP 2 > e > 0. 

In order to show that the propositions, say Pj and Pjk, are complementary, we compute 

[Pj>Pik] = h(\j)(k\-\k)V\)=iC (19) 

and (C) = ykj. Finally, 



\/pj( 1 -Pj)\/pjk(~L-Pjk) > \\Vkj\- (20) 

The variable ykj measures complementarity of Pj and Pjk- 

In supersymmetric projector model we start with 7 : 7i —>■ 7i and define the 'supercharge' 
Q :H®H^H®H, 

n ( ° A r* (n* \ Q 2 

^ 7 t oy \ 7 t 7 y Tr (Q 2 ) 

The situation where there are two types of super-partner density matrices is encountered in 
semantic analysis j^J. In semantic contexts the two density matrices correspond to words 
and sentences, respectively, and it is an intriguing issue if a link to DNA is here accidental. 
We use the model to describe the process of unwinding of two helices. 



Any of the probabilities we have discussed above is consistent with the frequency inter- 
pretation of experiments since always < p < 1. The subleties related to different models 
become visible when one comes to joint and conditional probabilities, since only with at 
least two propositions the issues of non-commutativity become relevant. 

The act of conditioning by a positive answer to the question represented by a projector 
P changes density matrices according to 

>" (22) 

Joint probability of successes in a sequence of measurements P^ P 2 , . . . , P k reads 

p(Pi n ■ ■ • n P k ) = Tr(p 1 ...P k -. 1 p k p k _ 1 ...p lP ). 

Here Pi is the first measurement and P k — the last. Since projectors belonging to differ- 
ent maximal sets do not all commute with one another, the joint probabilities have to be 
associated with noncommutativity of some sort. In molecular context the most natural non- 
commutativity is associated with chemical reactions: It is typical that the reactions X — > Y 
and Y — > X involve different kinetic constants. In particular, it is possible that the chain 
X — > Y — > Z of elementary reactions is possible, whereas the direct reaction X — > Z is not 
allowed. One can see here an analogy to the well known Malusian experiment with photons 
and three polarizers. It is worth mentioning that the very fact that some diagrams of the 
form X — > Y are not allowed may imply a Hilbert space structure j^j] , of course under some 
asumptions about the dynamics. 



III. SOLITON KINETIC EQUATIONS 



Let us characterize soliton kinetic equations in a general way. In explicit examples we 
will work only with a particular case, but since the subject is not widely known we find it 
important to set our results in a more general context. For further technicalities the readers 



are referred to 



22, 



231 ] . Consider an operator- valued function 



M 



X X = J2 A n A (n) + (23) 

n=0 n=l 

of a complex spectral parameter A. An invertible operator D x is a Darboux matrix for X\ if 

L M 

X X [1] = D~'X X D X = ^ X n X (n) [1] + Y, ^ n X { - n) [1] (24) 



n=0 



n=l 



where the X( n )[l], X(_ n )[l] are A-independent 
with 



24, 



26 



27 



28, 



Dy 



v — a 

i + — v p 

jj, — A 



. In this paper we work 



(25) 



where P = P 2 is a projector constructed in terms of right and left eigenvectors of, respec- 
tively, and X v [20, |3jJ, If v = fj, and P = then D\ is unitary: D A 1 = Dy 
Of particular importance is the case 

L 
n=0 

since then putting A = in (J2fijl we obtain 

X (0) [l] = X Q [1] = D^XqDq = D Q l X m D Q . 



(26) 



(27) 



In later applications to kinetic equations we will work with unitary D\ and Xq = p. 
Now consider two operator-valued functions F(X) and G(X) restricted only by 

[F(X),X] = [G(X),X] = 0, 
FiTXT- 1 ) = TF(X)T-\ 
G(TXT~ l ) = TG(X)T-\ 



for any T. The Lax pair is 



z\ip\ = ^aA(a), -iipx = if>\Y(\) 



where Z\ are complex numbers, and 

N 

= ^A fe X fc , 

/c=0 
L 



L M 



k=0 



A, 



Bi 



(L-k)\ \ds L 



k=l 
d L ~ k 



A fe 



c=0 



d 



M-k 



- k G(X [e) ) 



6 = 



(M-k)\ \de M ~ 

The compatibility condition for the Lax pair consists oniV + 1 equations 

N m 

%X m = A m-k] + X)[*k, B k - m ], 0<m<N. 



(28) 
(29) 
(30) 

(31) 

(32) 

(33) 
(34) 
(35) 

(36) 



k=m+l 



k=0 



For k < the operators are defined by means of (JMj) an d ()35|). The equation for m = is 
the von Neumann equation in question, and those for m > are additional relations between 
other operators occuring in the von Neumann equation. 

The Darboux matrix defines a Darboux transformation for the Lax pair and its compat- 
ibility conditions as follows. Take A( n )[l] defined via (|2IJ), and A fc [l], B k [l] defined by (JHlj) . 
(jnU), but with X( n ) replaced by A( n )[l]. Then 

N m 

iX m [l) = Ml], + Y^[X k [l], B k - m [l}}, < m < A, (37) 

k=m+l k=0 

and 

*aV>a[1] = ^a[1]X (a) [1], -# a [1] = Va[1]V(a)[1] (38) 

with ^a[1] = ip\D\. Restricting the operators X m to matrices and defining for X = p = 
> the kinetic variables according to the recipes of Sec. II we arrive at systems of soliton 
kinetic equations. 

IV. ONE-FIELD EQUATIONS 

We now restrict the explicit examples to 'one-field' equations, that is, those corresponding 
to iV = 1, i.e. 

X w = X + \X 1 = p + \H. (39) 

The equations are one-field since only X = p depends on time. One of the two dynamical 
equations reads here X\ = 0. Examples of N > 1 equations can be found in 2, An 
interesing class of nonlinearities occurs for F(X) — 0, M — 1, G(X) = —f(X). Then 

ip=[H,f(j>)]. (40) 

The right-hand-side of (J4*U|) can be rewritten as [Hf(p),p], with Hf(p) whose explicit form 
depends on /, showing that this is indeed a von Neumann equation. The equation itself 
was introduced in j^] and later discussed in the context of non-extensive statistics in 3(|. 
The soliton technique of integration based on Darboux-covariant Lax representation was 
introduced in and a class of 1-soliton self- switching solutions was there found for f(p) = 
p 2 . The Lax-pair representation and Darboux covariance properties of the general form were 
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established in |38j and explicit self- switching solutions were there found for f(p) = p q —2p q 1 , 
where q is an arbitrary real number. The explicit form of the Lax pair is simple 



iipx = j4>\f(p), z x i)\ = V>a(p + A-H"). 



(41) 



In the context of autocatalytic networks the examples of interest involve third order polyno- 
mial nonlinearities: f(p) = a 3 p 3 + a 2 p 2 + a\p. Let us choose f(p) = p 3 — 2p 2 . The equation 
can be written in two equivalent forms 



ip = [H,p 3 -2p 2 ] = [H(p),p], 
H(p) = p 2 H + pHp + Hp 2 - 2{Hp + pH) 



(42) 
(43) 



Taking 



^0 ^ 



H 



k 
2k 



(44) 



and performing the splitting (JT5j) we obtain six coupled kinetic equations with nonlinearities 
of third order. A closer look at their structure reveals that the reduction of variables from 
six to four can be obtained if one imposes the constraint 



P23 = ZP12, P32 = Zp21, pn = P33, 



1. 



(45) 



Assuming the above reduction with z — — 1 and introducing the four variables, 



.4 
B 
C 
D 



[A] ~ 2 Pl - 2 P2 



[B] 

^ 2 
[D)-Pi 



Pi = aria 

1 1 

Pi ~ 

2/13, 



Xl2 



2/12 



(46) 
(47) 
(48) 
(49) 



where and are, respectively, real and imaginary parts of pij, [A], [B], [C], [D], are 
nonnegative for all times and p\ = pn, p 2 = P22 are constants, we can rewrite the von 
Neumann equation as the set of four kinetic equations 
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[A] = -k [2C 3 + 2A 2 C - AD(2 Pl + P2 -2) 

+ C(B 2 + D 2 - 2 Pl +pl- 2p 2 + PlP2 +p 2 2 + B(2 Pl + p 2 - 2))) 

[B] = -2k (L> 3 + 2A 2 D + D(B 2 + 2C 2 - A Pl + Spf) - 2AC(2 Pl +p 2 -2 

[C] = k(2A 3 -CD(2 Pl +p 2 -2) 
+ A(B 2 + 2C 2 + D 2 — 2 Pl + p\- 2 P2 + PlP2 + p 2 2 - B(2 Pl + P2 - 2 

[D] = 2k(B 3 + B(2C 2 + D 2 + Pl (-4 + 3 Pl )) 
+ A 2 (2 + 2B- 2 Pl - P2 ) + C 2 (2 Pl + P2 -2 



(50) 
(51) 

(52) 

(53) 

It is possible to derive this system of kinetic equations from a system of elementary reactions 
of order not higher than three. This system of equations is a particular case, for q = 3, of 

;he general equation ip = [H,p q — 2p q ~ 1 ]. Its solution for an arbitrary value of q, found in 

38| by means of Darboux transformations, reads 



m 
at) 



( 3/2 m ^ 
-m 7/4 at) 
\ at) at) 3/2 j 

(Si + Vfye"^-^ 
AV2(l + e 2 ^(*-*°)) C 
1 - iV3 - 2e 2w ^- to ^ 



ikt 



,2ikt 



(54) 

(55) 
(56) 



4(1 + e2u q (t~ta)\ 

oj q = [(4/7) 1-9 — l]k/y/3. This solution is not normalized and hence the kinetic variables 
have be to rescaled if one is interested in probablities: P a = [A]/ Trp, and so on. Properties 
of the solution are illustrated by the figures. The scale of simplification and the efficiency of 
the von Neumann formalism become evident if one compares the von Neumann and kinetic 
forms of the same dynamics. 



V. FLUCTUATIONS VIA DARBOUX TRANSFORMATIONS 

The set of solutions of von Neumannn kinetic equations is naturally equipped with the 
Hilbert-Schmidt norm. The distance between two solutions p(t) and p(t) is 

|| p(t) - pit) ||= V 'Tr (p(t) - p(t)) 2 '. (57) 
12 



The solutions (|54j) describe families of distorted helices. We will now explicitly show that for 
any distorted helix (|54 )) -([56 p there exist two families of undeformed helices corresponding 
to two different solutions p±(t) of the same equation, such that (i) eigenvalues of p(t), p+(t), 
and p-(t) are identical, (ii) p(t) and p±(t) are linked by Darboux transformations, (iii) pit) 
and p±(t) belong to the same symplectic leaf of the Lie-Poisson structure corresponding to 
the dynamics (i.e. their Poisson-reduced dynamics occurs in the same phase-space), and (iv) 

\\p(t)-p±(t) \\<e (58) 

for any e > and for an infinitely long period of time t. In other words, an arbitrarily 
long part of an undeformed helix is arbitrarily close to a helix that involves a defect. In 
principle, a small fluctuation of the undeformed helix, at some t and along a symplectic 
leaf of the dynamics, may map it into a helix with a defect. However, if one restricts the 
analysis to a finite-length helix t\ < t < £2 (corresponding to a finite number of steps 
t n = nAt), then the defect may or may not become visible, because different Darboux- 
transformation fluctuations of the same p(t) result in helices with different locations of the 
defect (parametrized by to)- How quickly will the defect occur depends on what kind of 
a Darboux transformation is responsible for the fluctuation p(t) 1— > p(t). To see how this 
works consider 



( 3/2 ( ± (t)\ 



P±{t) 



Ut) = -\e™\ C-{t) = l -^e™. (59) 



7/4 
\C±(t) 3/2 J 

Density matrices p±(t) are the asymptotic forms of p(t) for —u> q to — > ±00. These two q- 
independent asyptotes also satisfy ip± = [H,p± — 2p q ± l ] for any real q, which follows from 
the fact that they can be obtained from the same solution by Darboux transformations 
differing by initial conditions for the Lax pair. This immediately implies that all these 
density matrices are isospectral and possess the same values of the Casimirs Tr(p n ). The 
Casimirs allow to reduce the Lie-Poisson structure to the same symplectic leaf, and thus 
the fluctuations occur in the same phase-space. The corresponding Hamiltonian function is 
given by Trif/(p). 
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VI. FORMATION OF OPEN STATE IN THE SUPERSYMMETRIC MODEL 



The supersymmetric model naturally leads to the double-helix structure. Consider an 
arbitrary operator 7 whose relation to a density matrix is p = 7^7. A superpartner of p is 
defined in the usual way as p = 77^. This type of structue is encountered in Latent Semantic 
Analysis where p and p represent word- vector and sentence- vector density matrices . 

We know that the Darboux transformation acts on density matrices by p 1— > p\ = DqpD . 
This is equivalent to 7 1— > 71 = jD , if p\ = 7^71 , or the supercharge transformation 




7 f / \ D I \j{0 

where Q is given by (J21|) . The action of the Darboux transformation on the super-density 
matrix is 

\ 7i f 7i J \ Dy-yDo J \ Pl J 
since D is unitary. At the level of super-density matrices the Darboux transformation acts 
by 

/ n n \ / n n \ 

(61) 

Of particular interest is the symmetric case 77^ = 7^7. This is trivially satisfied if 7 = 7"^ 
and 7 itself is a supercharge for p = p. Such a 7 always exists since p, being a positive 
operator, possesses a unique square root yfp = 7. The Darboux transformation breakes this 
symmetry as follows 





(62) 



P 0\ P 

Op/ \ Pl 

It is interesting that any sequence of Darboux transformations is here necessarily cyclic in 
the sense that 

/ po\f p \ ^ ^ I p o\f P o 

p) ~* \0 Pl ) ~* ~*\0p n )^\0 P 

with n being related to the dimension of the matrix p. To understand why this has to 
happen let us recall that the Darboux transformation involves two solutions of a Lax pair: 
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The one that is transformed, and the one that is used to define the transformation. The 
transformation is trivial if the two solutions are linearly dependent. It is clear that the 
number of linearly independent solutions is limited by the dimension of the linear space 
they live in. In the context of 3 x 3 solutions we find 

P -> Pi -> P2 -> Ps -> P, (63) 

no matter what (nontrivial) transformations have been performed, and two of these four 
density matrices have the asymptotic form p±(t). This is why after four Darboux fluctuations 
the system returns to its original state. 

For super-density matrices this process is strikingly similar to what happens with DNA 
during formation of an open state. Initially we have two identical strands described by the 
symmetric state p = p. The first nontrivial fluctuation maps it into the pair (p,pi) whose 
evolution looks as either Fig. 1 or Fig. 2 where p and p\ correspond to the two helices. Let 
us note, that the helices are practically overlapping over arbitrarily long parts, and therefore 
an arbitrarily small fluctuation may turn one into another. If the fluctuation is of the form 
shown in Fig. 1 then somewhere in the future the two helices will unwind, unless a new 
fluctuation will prevent it. If the fluctuation is of the form shown in Fig. 2, then the helix 
switches into a helix that wants to return to the asymptote, and we will not see a change of 
form. However, if the helices disentangled, then the helix approaches the other asymptote, 
and a subsequent small fluctuation will not be observed. Finally, another fluctuation maps 
it into a state that starts to return asymptotically to the original state. 

The fluctuations deform the helices according to the scheme (|63j) . However, we do not 
know which Darboux transformation is responsible for p — > pi, say. The point is that this 
is not really important! Fig. 5 shows two helices corresponding to two randomly chosen 
Darboux transformations. The random element is when do the helices start to unwind. 
All these helices have the same asymptotic forms and the part where the two helices look 
different is the one where the Hilbert-Schmidt distance is maximal. 

One can illustrate this phenomenon if one plots the Hilbert-Schmidt distance between 
two helices whose forms differ by the parameter to- Fig. 6 shows the plot of the Hilbert- 
Schmidt distance between two solutions that were produced from the same seed solution 
by two different Darboux transformations, as a function of t and the fluctuation parameter 
to (in one solution to = and is arbitrary in the other one). Two such density matrices 
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effectively differ from each other in a finite time interval. 

VII. SUMMARY AND PERSPECTIVES 

We have started with the observation that a system evolving into helical structures has to 
resemble an oscillator, at least at a sufficiently abstract level of description. Armed with this 
intuition and the fact that dynamical equations of DNA must be kinetic in some sense we 
have arrived at equations of a von Neumann type. Then, imposing soliton-type integrability, 
we restricted the class to Darboux-covariant von Neumann equations. 

The formalism we have proposed has many possible continuations. We have touched only 
the tip of the iceberg. Even the one-field equations are more general than those with the 
right-hand side [H, f(p)]. As shown in j23j], the class contains also Toda and other integrable 
lattices as particular cases. Solutions can be found by the same Darboux transformation, but 
one has to employ a different explicit technique of finding nontrivial seed solutions (cf. |4fl|). 
Still another technique of finding appropriate seed solutions that generalizes to arbitrary 
H and f(p) the method we have worked with, was introduced in [41| . When applied to 
explicit matrix examples the method leads again to qualitatively the same solutions as those 
we already know. One can risk the statement that the type of behavior (unwinding and 
switching between a finite number of states) we have described is generic for the whole class 
of integrable kinetic equations. In the context of genericity one should mention also the 
stability of solutions for f(p) = p q — 2p q ~ 1 under modifications of q (cf. the discussion in 
3^|). This means that the behaviour of the helix is basically the same for q = 3 and, say, 
q = 1000, or a non-integer q. If one tried to write down kinetic equations corresponding 
to q — 1000 one would have found the set of equations with polynomials of order 1000 and 
thousands of terms. Analogously, there would have been thousands of elementary very-high 
order reactions that, effectively lead to almost identical replication and unwinding as q = 3 
or q = 2. This fact may explain why incredibely simplified kinetic schemes one employs 
in replication retworks lead to evolutions that reasonably approximate various effects one 
expects in DNA-type systems. 
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FIG. 1: Formation of open state. The red, black, blue, and green helices represent, respectively, 
the matrix elements p-(t)is, /o(i)i3, /0_(i)i2, and p{t)\2- p(t) and p-(t) are given by and 
()59|) , The helices practically overlap in the past and unwind in a neighbourhood of t = 0. The 
four helices behave as if the external ones were the tapes, and the internal ones represented the 
read/write head. 
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